#include <iostream>
#include "na.h"

using namespace std;

int main()
{
 Tna ivp;
 *((Tpoint_kin_params*)&ivp)=TR;
 double t,te=10;
 ivp.dt=0.1;
 ivp.ro=0.1*ivp.betta;
 ivp.n1=ivp.n0=1;
 ivp.x0();
 for(t=0;t<te;t+=ivp.dt)
 {
  ivp.step();
  cout<<t<<' '<<ivp.n1<<endl;
 }
 return 0;
}